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We present measurements of the semileptonic decays B~ —* D°t~V t , B~ — ► D*°t~V t , B° — > 
D + r~V T , and B° — > D* + r~T7 T , which are potentially sensitive to non-Standard Model amplitudes. 
The data sample comprises 232 x 10 6 T(4S) — ► BB decays collected with the BABAR detector. 
From a combined fit to B~ and B° channels, we obtain the branching fractions B(B — > _Dt~77 t ) = 
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(0.86_± 0.24 ± 0.11 ± 0.06)% and B(B -> D*t~V t ) = (1.62 ± 0.31 ± 0.10 ± 0.05)% (normalized for 
the B°), where the uncertainties are statistical, systematic, and normalization-mode-related. 

PACS numbers: 12.15.Hh, 13.20.-v, 13.20.He, 14.40.Nd, 14.80.Cp 



Semileptonic decays of B mesons to the r lepton — 
the heaviest of the three charged leptons — provide a 
new source of information on Standard Model (SM) pro- 
cesses [1M2J3L as weu as a new window on physics beyond 
the SM^JalStzLli]. In the SM, semileptonic decays oc- 
cur at tree level and are mediated by the W boson, but 
the large mass of the r lepton provides sensitivity to ad- 
ditional amplitudes, such as those mediated by a charged 
Higgs boson. Experimentally, b — > cr~V T decays [§| are 
challenging because the final state contains not just one, 
but two or three neutrinos as a result of the r decay. 

Branching fractions for semileptonic B decays to r 
leptons are predicted to be smaller than those for I — 
e, \x [10j. Calculations based on the SM predict B(B° — > 
D+t-V t ) = (0.69 ± 0.04)% and B{B Q -> D*+ T -V T ) = 
(1.41±0.07)% [H, which account for most of the predicted 
inclusive rate B{B -> X cT -V T ) = (2.30±0.25)% 12] (here, 
X c represents all hadronic final states from the b — » c 
transition). Calculations 0, d, @, 0, Hi m multi-Higgs 
doublet models show that substantial departures, either 
positive or negative, from the SM decay rate could oc- 
cur for B(B Dt~V t ). Those for B{B D* T -V T ), 
however, are expected to be smaller. 

Theoretical predictions for semileptonic decays to ex- 
clusive final states require knowledge of the form factors, 
which parametrize the hadronic current as functions of 
q 2 = (pb —Pdm) 2 - F° r hght leptons (e, /i), there is ef- 
fectively one form factor for B — > D£~Ve, while there 
are three for B — > D*£^V(. If a r lepton is produced 
instead, one additional form factor enters in each mode. 
The form factors for B — ► D^l~Th decays involving the 
light leptons have been measured [Tl| . Heavy quark sym- 
metry (HQS) relations 12| allow one to express the two 
additional form factors for B — ► D(*'t~V t in terms of 
the form factors measurable from decays with the light 
leptons. With sufficient data, one could probe the addi- 
tional form factors and test the HQS relations. 

The first measurements of semileptonic 6-hadron de- 
cays to t leptons were performed by the LEP exper- 
iments [l3| operating at the Z° resonance, yielding 
an average [14| inclusive branching fraction B(bhad — > 
Xt~v t ) = (2.48±0.26)%, where 6had represents the mix- 
ture of 6-hadrons produced in Z° — ► bb decays. The Belle 
experiment has recently obtained B(B° — > D* + t~V t ) = 



(2.02±°;« ± 0.37)% [l 

Wc determine the branching fractions of four exclu- 
sive decay modes: B~ — ► D°t~T7 t , B~ — > D* t~17 t , 



B u 



D+t~V t , and B° 



D 



* + t~T> 



each of which is 



measured relative to the corresponding e and /1 modes. 
To reconstruct the r, we use the decays r~ — * e~v e v T 
and r~ — > yTv^Vr, which arc experimentally most ac- 



cessible. The main challenge of the measurement is to 
separate B — ► D^*>t~V t decays, which have three neu- 
trinos, from B — ► D^*>l~Vi decays, which have the same 
observable final-state particles but only one neutrino. 

We analyze data collected with the BABAR detector [l6[ 
at the PEP-II e+e~ storage rings at the Stanford Lin- 
car Accelerator Center. The data sample used com- 
prises 208.9 fb _1 of integrated luminosity recorded on 
the T(4S) resonance, yielding 232 x 10 6 BB decays. 

The analysis strategy is to reconstruct the decays of 
both B mesons in the T(AS) — > BB event, providing 
powerful constraints on unobserved particles. One B 
meson, denoted Stag, is fully reconstructed in a purely 
hadronic decay chain. The remaining charged particles 
and photons are required to be consistent with the prod- 
ucts of a b — > c semileptonic B decay: a hadronic sys- 
tem, a meson, and a lepton (e or fi). The lep- 
ton may be either primary or from t~ — > l~V(V r . We 
calculate the missing four-momentum p nllss — [p e + e - — 
Ptag — Pd(*> ~ Pi] °f an y particles recoiling against the 
observed Stag + D^*'t system. A large peak at zero in 
m miss = Pmiss corresponds to semileptonic decays with 
one neutrino, whereas signal events form a broad tail out 
to m^ iss ~ 8 (GeV/c 2 ) 2 . To separate signal and back- 
ground events, we perform a fit to the joint distribution 
of m^jgg and the lepton momentum (|p||) in the rest 
frame of the B meson. In signal events, the observed 
lepton is the daughter of the r and typically has a soft 
spectrum; for most background events, this lepton typi- 
cally has higher momentum. 

We reconstruct B tag candidates [13] in 1114 final states 
-Stag — * D^Y ± . Tag-side candidates are recon- 

structed in 21 decay chains, and the Y system may 
consist of up to six light hadrons (7r ± , it , K ± , or Kg). 
Stag candidates are identified using two kinematic vari- 
ables, m E s = \A/4- |pta g | 2 and AE = E tag - y/s/2, 
where is the total e + e~ energy, |ptag| is the magni- 
tude of the Btag momentum, and Etag is the Bt ag energy, 
all defined in the e + e~ center-of-mass frame. We require 
TO ES > 5.27 GeV/c 2 and \AE\ < 72 MeV, correspond- 
ing to ±4cr (standard deviations). We reconstruct B ta , g 
candidates in approximately 0.3% to 0.5% of BB events. 

For the B meson decaying semileptonically, we re- 
construct £)W candidates in the modes D° — > K~tt + , 
K-ir+n , K-iT+n+ir-, K a s ir+ir~; D+ -> K--K+ir+, 
K-tt+tt+tt , K°ir + , K~K + -k + ; D*° -> D°ir°, D° T , 
and D*+ -> D 6 tt+, D+tt°. D (£>*) candidates are se- 
lected within 4(7 of the D mass (D* ~ D mass differ- 
ence), with cr typically 5-10 MeV/c 2 (1-2 MeV/c 2 ). Elec- 
tron candidates must have lab- frame momentum |p e | > 
300 MeV/c; muon candidates must have an appropriate 
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signature in the muon detector system, effectively requir- 
ing |p M | ~ 600MeV/c. The energy of electron candidates 
is corrected for bremsstrahlung energy loss if photons are 
found close to the electron direction. 

We require that all charged tracks be associated with 
either the B ts , g , D", or I candidate. We compute -Eextra, 
the sum of the energies of all photon candidates not as- 
sociated with the £> tag + D^*>£ candidate system, and we 
require -E ex t ra < 150-300 MeV, depending on the DW 
channel. We suppress hadronic events and combinato- 
rial backgrounds by requiring |p m i ss | > 200 MeV/c and 
q 2 > 4 (GeV/c 2 ) 2 . If multiple candidate systems pass 
this selection, we select the one with the lowest value of 
-Eextra- To improve the TO 2 ^ resolution, we perform a 
kinematic fit to the event, constraining particle masses 
to known values and requiring tracks from B, D, and 
K® mesons to originate from appropriate common ver- 
tices. All event selection requirements and fit procedures 
have been defined using simulated events or using control 
samples in data that exclude the signal region. 

Figure [T] shows the distributions of m^iaa f° r the four 
D^*'£ channels, along with the projections of the max- 
imum likelihood fit discussed below. We observe large 
peaks at TO 2 ^ « as well as events in the signal region at 
large m^^. The peaks are mainly due to B — > D^*'l~Vi, 
which serve as normalization modes. The structure of 
this background is shown in the inset figures, which ex- 
pand the region -0.4 < < 1.4 (GeV/c 2 ) 2 . B — > 
D*£~Ve background is the dominant feature in the two 
D*£ channels (Figs. UK, c); the two D£ channels (Figs.QJ), 
d) are dominated by B — > D£~Vg decays but also include 
substantial contributions from true D* mesons where the 
low-momentum 7r° or photon from D* — > Dir° or is 
not reconstructed. Similarly, B — > D*t~T> t events can 
feed down to the D£ channels. The fit therefore includes 
feed-down components for both the signal and normal- 
ization modes, as well as smaller feed-up contributions 
from B — ► D{£~~ /t~)V into the D*i channels. Other 
sources of background include B — > D**(£~ /t~)v events 
(here D** represents charm resonances heavier than the 
D*(2010), as well as non-resonant D^rnr systems with 
n > 1); charge-crossfeed (B — > D^£~T>i events recon- 
structed with the wrong charge for the J5 ta g and D(*> me- 
son, typically because a low-momentum is swapped 
between them); and combinatorial background. This last 
background is dominated by hadronic B decays, such as 
B —>■ D^Di*\ that produce a secondary lepton, includ- 
ing t leptons from D s decay. 

To constrain B — ► D**(£~ /t~)v background, we select 
four control samples, identical to the signal channels but 
in which an extra tt° meson is observed. Most of the D** 
background in the signal channels occurs when the 7r° 
from D** — ► D^t: is not reconstructed, so these con- 
trol samples provide a normalization of the background 
source. D** decays in which a tt^ is lost do not have the 
correct charge correlation between the B tag and D", 




mLs(GeV/c 2 ) 2 

FIG. 1: (Color online) Distributions of events and fit projec- 
tions in m^jjgg for the four final states: D*°£~ , D°£~ , 
and D + £~ . The normalization region m^ iss w is shown with 
finer binning in the insets. The fit components are combina- 
torial background (white, below dashed line), charge cross- 
feed background (white, above dashed line), the B — > Dl~Ve 
normalization mode (// hatching, yellow), the B — > D*l~T>t 
normalization mode (\\ hatching, light blue), B — > D**l~Vi 
background (dark, or blue), the B — > Dt~v t signal (light 
grey, green), and the B — > D*t~V t signal (medium grey, ma- 
genta). The fit shown incorporates the B~-B° constraints. 

and decays with two missing charged pions are rare. The 
feed-down probabilities for the D**{£~~ /t~)V background 
are determined from simulation, with uncertainties in the 
D** content treated as a systematic error. 

We perform a relative measurement, extracting both 
signal B — > D(*^t~V t and normalization B — > D^*'£~Vg 
yields from the fit to obtain the four branching ratios 
R(D°), R(D+), R(D*°), and R(D*+), where, for exam- 
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pie, R(D*°) = B(B~ -> D*°t-V t )/B(B- -> D^fVi). 
These ratios are normalized such that £ represents only 
one of e or fi; however, both light lepton species are 
included in the measurement. Signal and background 
yields are extracted using an extended, unbinncd max- 
imum likelihood fit to the joint (m^ iss , |pj|) distribu- 
tion. The 18-parameter fit is performed simultaneously 
in the four signal channels and the four D** control sam- 
ples. In each of the four signal channels, we describe the 
data as the sum of seven components (shown in Fig.[T]): 
Dt-V t , D*t~V t , Di~V t , D*£~V e , D**(£-/t~)V, charge 
crossfeed, and combinatorial background. The four D** 
control samples are described as the sum of five compo- 
nents: D**(£~ /t~)V, D£~17e, D*£~~p£, charge crossfeed, 
and combinatorial background. Probability distribution 
functions (PDFs) are primarily determined from simu- 
lated event samples. Both the signal and normalization 
modes are described using HQET-based form factors [HI 
for which the parameters and their uncertainties are de- 
termined by experimental measurements 11]. Param- 



modes contribute 0.3% or less. Finally, the uncertainty 
on B(t~ — > £~~Vf.v T ) [l_4| contributes 0.2% to all modes. 



eters describing the amount of the dominant feed-down 
components — D* feed-down into the D channels — are de- 
termined directly by the fit. 

Table H] summarizes the results from two fits, one in 
which all four signal yields can vary independently, and 
a second fit in which we constrain [19( R(D + ) — R(D°) 
and R(D* + ) = R(D*°). The m^ iss projections shown in 
Fig. [T]are those from this i? _ -B°-constrained fit. 

The features of the event sample have been exten- 
sively checked. The observed lepton spectra are well 
described by the fit both in signal- and in background- 
dominated regions. The properties of reconstructed £? tag 
mesons, such as charged and neutral daughter multiplic- 
ities, are consistent with expectations. Control samples 
of B — > D^£~V( events, kinematically selected without 
a cut on m^jgg, provide checks of numerous distributions, 
including the TO^ iss tails. 

Systematic uncertainties on R associated with the fit, 
(AR/ R)nt in Table[lJ are determined by running ensem- 
bles of fits in which input parameters are distributed ac- 
cording to our knowledge of the underlying source, and 
include the PDF parametrization (2% to 12%); the com- 
position of combinatorial backgrounds (2% to 11%); the 
mixture of D** states in B — > D**£~v t decays (0.3% to 
6%); the B -> D* form factors (0.2% to 1.9%); the vr° 
reconstruction efficiency, which affects the D* and D** 
feed-down rates (0.5% to 1.1%); and the m;^ iss resolution 
for B -> D*£~Vi events (0.1% to 1.6%). Uncertainties 
on the B — > D form factors contribute less than 1%. Un- 
certainties on R propagated from the ratio of efficiencies, 
(AR/R) E in Table [H are typically small due to cancel- 
lations, and include the limited statistics in the simula- 
tion (0.8% to 1.5%) and systematic errors related to de- 
tector performance (0.2% to 0.7%). Uncertainties from 
modeling final-state radiation are 0.3% to 0.5%; uncer- 
tainties on the branching fractions of the reconstructed 



Table U gives the significances of the signal yields. The 
statistical significance is determined from ^/2A(ln£), 
where A(ln£) is the change in log-likelihood between 
the nominal fit and the no-signal hypothesis. The to- 
tal significance is determined by including (Ai?/i?)g t in 
quadrature with the statistical error. 

We have presented measurements of the decays B — > 
Dt~V t and B — > D*t~T> Ti relative to the corresponding 
decays to light leptons. We find R(D) = (41.6 ± 11.7 ± 
5.2)% and R(D*) = (29.7 ± 5.6 ± 1.8)%, where the first 
error is statistical and the second is systematic. Normal- 
izing to known B° branching fractions [23], we obtain 



B(B -> Dt-V t ) 
B(B -> D*t-V t ) 



(0.86 ±0.24 ±0.11 ±0.06)% 
(1.62 ±0.31 ±0.10 ±0.05)%, 



where the third error is from that on the normaliza- 
tion mode branching fraction. The significances of the 
signals are 3.6cr and 6.2<r, respectively. The modes 
B~ -► D°t-V t , B~ -> D*°t-V Ti and B Q -► D + T -V T 
have not been studied previously, while the measurement 
of B — > D* + t~V t is consistent with the Belle result [15]. 
The averaged branching fractions are about la higher 
than the SM predictions but, given the uncertainties, 
there is still room for a sizeable non-SM contribution. 
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D*t~V t 
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